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Abstract 



A Ricmannian metric with a local contraction property can be used to prove 
existence and uniqueness of a periodic orbit and determine a subset of its basin 
of attraction. While the existence of such a contraction metric is equivalent to 
the existence of an exponentially stable periodic orbit, the explicit construction 
of the metric is a difficult problem. 

In this paper, the construction of such a contraction metric is achieved 
by formulating it as an equivalent problem, namely a feasibility problem in 
semidefinite optimization. The contraction metric, a matrix-valued function, is 
constructed as a continuous piecewise affine (CPA) function, which is affinc on 
each simplex of a triangulation of the phase space. The contraction conditions 
are formulated as conditions on the values at the vertices. 

The paper states a semidefinite optimization problem. We prove on the one 
hand that a feasible solution of the optimization problem determines a CPA 
contraction metric and on the other hand that the optimization problem is 
always feasible if the system has an exponentially stable periodic orbit and the 
triangulation is fine enough. An objective function can be used to obtain a 
bound on the largest Floquet exponent of the periodic orbit. 

1 Introduction 

In this paper we consider a time-periodic ODE of the form x = f(t,x), where 
f{t, x) = f(t + T, x) for all (t, x) G R x R n with a fixed period T > 0, and study the 
basin of attraction of a periodic solution. 
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The basin of attraction can be computed using a variety of methods: Invariant 
manifolds form the boundaries of basins of attraction, and their computation can 
thus be used to find a basin of attraction P33 Hj- However, this method needs 
additional arguments to ensure that a certain region is the basin of attraction of an 
attractor, and that, for example, there are no other attractors in that region. Other 
approaches to compute the basin of attraction are for example the cell mapping 
approach [T7] or set oriented methods [6j which divide the phase space into cells 
and compute the dynamics between these cells. 

Lyapunov functions [21] are a natural way of analysing the basin of attrac- 
tion, since they start from the attractive solution, not from the boundary. Moreover, 
through their level sets, they give additional information about the basin of attrac- 
tion than just the boundary. Converse theorems which guarantee the existence 
of a Lyapunov function under certain conditions have been given by many authors, 
for an overview see |15j . However, all converse theorems offer no general method to 
analytically construct Lyapunov functions. 

Recently, several methods to construct Lyapunov functions have been proposed: 
Hafstein [TJ] constructed a piecewise affine Lyapunov function using linear pro- 
gramming. Parrilo [21] and Papachristodoulou and Prajna in [23] consider the nu- 
merical construction of Lyapunov functions that are presentable as sum of squares 
of polynomials (SOS) for autonomous polynomial systems. These ideas have been 
taken further by recent publications of Peet [25J and Peet and Papachristodoulou 
|26j , where the existence of a polynomial Lyapunov function on bounded regions for 
exponentially stable systems in proven. 

A different method deals with Zubov's equation and computes a solution of 
this partial differential equation (PDE). In Camilli et al. [SJ, Zubov's method was 
extended to control problems. Giesl considered a particular Lyapunov function 
satisfying a linear PDE which was solved using meshless collocation, in particular 
Radial Basis Functions [9] . This method has been extended to time-periodic ODEs 

Lyapunov functions attain their minimum on the attractor and have a negative 
orbital derivative for all points in the basin of attraction apart from the attractor. 
Hence, it is necessary to have exact information about the attractor in the phase 
space before one can compute a Lyapunov function. Whereas this information might 
be easy to obtain in special examples, in general this information is not available. 

Local contraction property — Borg's criterion 

Another method to characterise the basin of attraction, introduced by Borg 
[3], uses a local contraction property and does not require information about the 
periodic orbit. Let M(t,x) be a Riemannian metric, i.e. M G C^M x M. n , W ixn ) 
such that M(i, x) is a positive definite, symmetric (n x n) matrix for all (t, x). Then 
(v,w)m(x) '■= v T M{x)w defines a point-dependent scalar product, where v,w G M n . 



The sign of the real-valued function LM(t,x), cf. (1.1), then describes whether the 



solution through (t, x) and adjacent solutions approach each other with respect to 
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the Riemannian metric M. We define 



L,M(t,x) := max w 

ioeM n ,w T M{t,x)w=l 



T 



M(t,x)D x f(t,x) + ^M'(t,x) 



w, (1.1) 



where M'(t, x) denotes the orbital derivative of M(t, x), which is the derivative along 
solutions of x = f(t,x). 

If L]\f(t,x) is negative for all (t, x) 6 K where K is a positively invariant, con- 
nected set, then K is a subset of the basin of attraction of a unique periodic orbit 
in K. 



The maximum in (1.1) is taken over all w £ IR n wit h a norm condition, 



and LM(t,x) < is equivalent to LM(t,x) < given by (1.2), where Xmax(-) 
denotes the maximal eigenvalue of a symmetric matrix. Here, we use that 
w T M(t,x)D x f(t,x)w = w T D x f(t,x) T M(t,x)w. 

L M (t,x) := max w T \M(t, x)D x f(t, x) + D x f(t, x) T M(t, x) + M'(t, x)] w 

weiR n ,||w||=i 

= Xmax(M(t,x)D x f(t,x) + D x f(t,x) T M(t,x) + M'(t,x)) . (1.2) 

We seek to find a matrix- valued function M satisfying Lm(J:,x) < 0. This is 
equivalent to the condition that the symmetric matrix 

M(t, x)D x f(t, x) + D x f(t, x) T M{t, x) + M'(t, x) 

is negative definite. As this is a Linear Matrix Inequality, it can be formulated as a 
constraint of a semidefinite optimization problem. 

While the sufficiency of this local contraction criterion in the autonomous case 
goes back to [31 [161 123 120] , its necessity was shown in [7] . The method was extended 
to time-periodic systems [8]. 

The advantage of this method over, for example, Lyapunov functions, is that it 
does not require information about the position of the periodic orbit. Moreover, the 
criterion is robust to small errors. 

Although the existence of Riemannian metrics has been shown [7] , it remains a 
difficult problem to construct them for concrete examples. This is a similar problem 
to the construction of a (scalar-valued) Lyapunov function, but Borg's criterion re- 
quires the construction of a matrix- valued function M(t,x). In the two-dimensional 
autonomous case, however, there exists a special Riemannian metric of the form 
M(x) = e 2W ^I, where W is a scalar- valued function [7J. This can be used to find 
an approximation using Radial Basis Functions [10]. In higher dimensions, however, 
the existence of such a special Riemannian metric is not true in general [7J . In [11] , a 
combination of a Riemannian metric locally near the periodic orbit with a Lyapunov 
function further away was used, and the construction was again achieved by Radial 
Basis Functions. This method, however, heavily depends on information about the 
periodic orbit, which was obtained by a numerical approximation of the periodic 
orbit and its variational equation. 
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In this paper, we will develop a new method to construct a Riemannian metric to 
fulfill Borg's criterion, which will use semidefinite optimization and does not require 
any information about the periodic orbit. 

Semidefinite Optimization 

A semidefinite optimization problem for the variables y\ , . . . , y m is of the form 

minimize YJlLi °iVi 

subject to YT=i F iVi ~ F o = X h 0, 

where F{ are symmetric (N x N) matrices and X y means that the matrix X is 
positive semidefinite. 

The goal of this paper is to formulate the condition of a contraction metric as 
a semidefinite optimization (feasibility) problem. In a subsequent paper we will 
discuss the details of how to solve this problem efficiently. 

The main idea is to first triangulate the phase space. The Riemannian met- 
ric, i.e. the symmetric matrix M(t,x), will be expressed as a continuous piecewise 
affine (CPA) function, i.e. if M is given at the vertices (to,xo), . . . , (t n+ i, x n+ i) of 
a simplex, then M(t, x) = Y^i=o AjM(ij, Xi), where (t, x) = Y^i=o ^i{^i^ x i)- The 
conditions of Borg's criterion will become the constraints of a semidefinite optimiza- 
tion problem on the vertices of the triangulation, which will ensure the contraction 
property for all points in the simplices. 

In [I], a contraction metric is also constructed using semidefinite optimization. 
There are, however, three main differences to our approach: firstly, adjacent trajec- 
tories in all directions are contracted, whereas in our case the contraction takes place 
in the n-dimensional subspace W 1 of W l+ , but not in the time-direction. Thus, in 
our case, the attractor is a periodic orbit, whereas in their case, it is an equilibrium 
point. Secondly, and more importantly, the above paper transforms the construc- 
tion problem to a Linear Matrix Inequality and solves this using a sum-of-squares 
approach. This approach is used to prove global stability, i.e. the basin of attrac- 
tion is the whole space. The contraction metric is a polynomial function, and the 
system considered is assumed to be polynomial, too. In this paper, we study sys- 
tems which are not necessarily polynomial nor globally stable, and we triangulate 
the phase space to obtain a large subset of the basin of attraction. Lastly, we are 
able to prove that the semidefinite optimization problem is feasible if and only if the 
dynamical system has an exponentially stable periodic orbit, whereas in their paper 
the equivalence does not hold since the sum-of-squares condition is more restrictive 
than positive definiteness of matrices. 

The paper is organised as follows: In Section [2] we start with preliminaries, and 
in Section [3] we generalise the existing theorem for a smooth Riemannian metric M 
to a CPA (continuous piecewise affine) Riemannian contraction metric. We show 
that the existence of such a metric is sufficient to prove existence and uniqueness of 
a periodic orbit and to determine a subset of its basin of attraction. In Section [4j we 
describe the triangulation of the phase space into a simplicial complex and state the 
semidefinite optimization problem. Furthermore, we show that the feasibility of the 
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semidefinite optimization problem provides us with a CPA contraction metric. We 
also discuss a possible objective function to obtain a bound on the largest Floquet 
exponent. In Section [5j finally, we show that the semidefinite problem is feasible, if 
the dynamical system has an exponentially stable periodic orbit. 

Altogether, this paper shows that the problem of finding a contraction metric is 
equivalent to the feasibility of a semidefinite optimization problem. 



2 Preliminaries 

In this paper we consider a time-periodic ODE of the form 

x = f(t,x), (2.1) 
where f(t, x) = f(t + T, x) for all (t, x) G R x R n with a fixed period T > 0. We 
denote f(t,x) = f f(tx) ) an< ^ ^ = &i x )- We study the equation on the cylinder 

Sj< x R n as the phase space, where Sj> denotes the circle of circumference T. We 
assume that / G C 2 (S\ x E n ,R n ). If even / G C 3 (S^ x R n ,R n ) holds, then we 
obtain improved results, in particular higher order approximations; this is, however, 
not necessary to derive the main result. We denote the (unique) solution of the ODE 
with initial value x(to) = xo by x[t) =: Sf(to,xo) and denote (i + ip mod T,x(t)) =: 
St(to,xo) G Sj, x R n . Furthermore, we assume that the solution exists for all t > 0. 

We use the usual notations for the vector and matrix norms, in particular we 
denote by ||M|| 

max '■ — maxj ^'=1^ 7i the maximal entry of a matrix. S n denotes 
the set of all symmetric real- valued nxn matrices. For a symmetric matrix M G S n , 
^max(M) denotes its maximal eigenvalue, we write M >z if and only if M is positive 
semidefinite and M ^ if and only if M is negative semidefinite. The convex hull 
is defined by 

(n+l n+1 ^| 

co(x ,xi, . . . ,x n+1 ) := < ^AjXj: Xj > 0, ^ \j = 1 > . 



3 A CPA contraction metric is sufficient for a periodic 
orbit 

It was shown in [8j that a smooth contraction metric implies the existence and 
uniqueness of a periodic orbit and gives information about its basin of attraction. 
In this paper, we will seek to construct a CPA contraction metric, which is not of 
the smoothness required in the above paper. 

In this section we show that we can relax the conditions on the smoothness of M 
to cover the case of a CPA contraction metric. We will require that M is continuous, 
Lipschitz continuous with respect to x, and the forward orbital derivative exists. We 
will later relate this to the construction of the CPA contraction metric on a suitable 
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triangulation, and we will prove that such a CPA metric will satisfy all assumption 



that we make in this section (Lemma 4.7). The proof of the main theorem, Theorem 
3.4, will closely follow [S]. 

First we define a weaker notion of a Riemannian metric, which does not assume 
that M is differentiable, but only that the orbital derivative in forward time exists. 

Definition 3.1 (Riemannian metric) M is called a Riemannian metric for 



( fg.jp , if M G C°(S^ x W n ,8 n ) where M(t,x) is positive definite. Moreover, 
we assume that M is locally Lip schitz- continuous with respect to x, i.e. for all 
(to, xq) G S)p x R n there exists a neighborhood U C x W 1 such that 

\\f(t,x)-f(t,y)\\<L\\x-y\\ 

holds for all (t, x), (t,y) G U. Furthermore, we assume that the forward orbital 
derivative M' + (t,x) is defined for all (t,x) G <Sy x W 1 , where M' + (t,x) denotes the 
matrix 

M>(t,x) = lim M(Seit, X ))-Mjt, X ) 

We have the following statements for functions with a right-hand side derivative. 

Lemma 3.2 Let gi,g 2 G C°(M,M) be RHS differentiable, i.e. g'+(x) := 
lim /l _ s>0 + s( x + h }-9( x ) ex { s f s _ L e t G: M ^-M be differentiable. 
Then g\(x) • g2(x) and G(g(x)) are RHS differentiable and 

{gi ■ g2)' + {x) = (gi)+(x) ■ g 2 (x) + gi(x) ■ {g2)' + (x) (3.1) 
(Gog)' + (x) = G'(g(x))-g' + (x). (3.2) 

Let g G C°(I,M), where I C K is open, be RHS differentiable. Let [xi,X2] C I, 
and let g'+{C) d£, exist and be finite. 
Then we have 

rx2 



C-X2 

I g'+(0 d£ = g(x 2 ) - g(xi) 

J XI 



Proof: The first two statements follow directly from the usual proofs, replacing the 
limit with the RHS limit. The last statement is a result due to Lebesgue, formulated 
originally for Dini derivatives. □ 

Lemma 3.3 Let V(t, x) be a function which is locally Lipschitz with respect to x 
and let lim e _ >0+ n(^o) + e/fa,*o))-V^o,*o) ^ 

Then the orbital derivative Vl_(to,xo) exists and is equal to 

rr/ /■ \ ,. V(S e {to,x )) -V(t ,x ) 
V + (to,x ) := hm 

V((t ,x ) + ef(t , x )) - V(t , x ) 

= hm . 

6^0+ 8 
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Proof: Recall that f(t, x) 



f(t,x) 



and x = (t, x). We use that 



V(Se(t ,x )) - V(t , x ) V(S e (t , x )) - V((t , x ) + 9(1, f(t , x ))) 



V((t , x ) + 9(1, /(t , s ))) - V(to, x ) 



Note that due to the Lipschitz continuity \V(t, y) — V(t, x)\ < L\\y — x\\ we have for 
small 9 



V(S e (t , xq)) - V((t , x ) + 9(1, f(t ,x ))) 



< L 



\Sg{to, X ) - Xq- 9f(t , Xq) 



By the Mean Value Theorem there are h\(9), . . . , h n (9) E [0, 1] such that 

d 

(Sg(t ,x ) - x )i = 9^ t (Sl i{e)d (to,x Q )) i = 9fi(S ht{6)e (t ,x )). 



Hence, 



lim 

0->O+ 



\(Sg(t ,X ) -X - 9f(t ,X ))i 



< lim \fi(S h j e)e (t ,x )) - fi(t ,x )\ 
= 



since the solution and fi are both continuous. Altogether, we thus have 
V(S e (t ,xo)) - V((t ,x ) - 9(1, /(t , x ))) 



lim 

6>->0+ 



0. 



□ 



Theorem 3.4 Consider the equation x = f(t,x), where x G W 1 , and assume that 
f G C°(5y x M n ,]R n ) and all partial derivatives of f order one with respect to x are 
continuous functions of (t,x). Let / G C x W 1 be a connected, compact and 
positively invariant set. Let M be a Riemannian metric in the sense of Definition 

roi 

Moreover, assume Ljif(t,x) < —v < for all (t,x) £ G, where 



LM(t,x) := sup w 

«i6t n ,io T M(i,x)iu=l 



M(t,x)D x f(t,x) + -M' + (t,x) 



w (3.3) 



Then there exists one and only one periodic orbit !1 C G which is exponentially 
asymptotically stable. Moreover, for its basin of attraction G C A(Q) holds. 
If 



J 

Jo 



p(t) T M' + (t,x(t))p(t) dt 
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exists and is finite for all solutions x(t) with x(0) G G and all functions p G 
C°([0, T],M n ), f/ien f/ie largest real part —uq of all Floquet exponents of Q satis- 
fies 

-u < -v. 

Note that in contrast to [8], here L,M(t,x) is not necessarily continuous, since 
M',(t,x) is not continuous in general. 

Proof: The only parts in the proof that need to be changed slightly from [8] are 
the proof of Proposition 3.1 and the estimate on the Floquet exponents; the other 
steps of the proof are exactly the same. 
Proposition 3.1 

We replace the temporal derivative in Step III of [H Proposition 3.1] by the 
forward temporal derivative of 

A{9) 2 := [S$(t ,x + V )-S%(t ,x )} T M(S e (t ,xo)) 

■[S$(to,x Q + ri)-S$(!to,x )]. 



Note that the product rule holds for the RHS derivative as usual, cf. Lemma 3.2 
Furthermore, we use the comparison lemma in the version for RHS limits; a more 
general version for Dini derivatives can be found in [18, Lemma 3.4]. This shows 
that the result of [HI Proposition 3.1] remains true for the Riemannian metric M as 
in Definitionl3.il 



Floquet exponent 

In this part of the proof in [8] we need to show that 

(In (p(t) T M(t,x(t))p(t)))' + dt 

= In (p(T) T M(T, x(T))p(T)) - In (p(0) T M(0, x(0))p(0)) 

where x(t) is the periodic orbit and p{t)e~ Uot with p(0) = p(T) ^ is a solution of 
the first variation equation y = D x f(t,x(t))y along the periodic orbit. Note that 
we have already used that the RHS derivative satisfies the product rule, cf. Lemma 



3.2 Using the same lemma, the composition with the differentiable function In is 
also RHS differentiable. To apply the last statement of this lemma, we need to show 
that 



/ {In (p(t) T M(t,x(t))p(t)))' dt 
Jo 



exists and is finite. Indeed, we have, using the chain rule in Lemma 3.2, that 

(p(t) T M(t,x(t)) P (t))' 



(In {p(t) T M(t,x(t))p(t)))' 



p(t) T M(t,x(t))p(t) 
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Since M(t, x(t)) is positive definite and continuous and p(t) 7^ 0, the denominator is 
bounded away from zero. We apply the product rule to the nominator and use that 

T 

p(t) T M' + (t,x(t))p(t)dt 

o 

exists and is finite by assumption, and that the other terms are smooth. This shows 
the theorem. □ 



4 A solution of the semidefinite optimization problem 
defines a contraction metric 

In this section we will state a semidefinite optimization problem and show that 
a feasible solution of the semidefinite optimization problem defines a contraction 
metric, satisfying the assumptions of Theorem |3.4| The solution will be a CPA 
matrix- valued function M, which is defined by the values at the vertices. 

We will first define a triangulation of the phase space and then define a Rie- 
mannian metric by its values at the vertices of this triangulation and affine on each 
simplex. This CPA Riemannian metric is shown to fulfill all conditions of Theorem 



3.4, also inside the simplices. 



4.1 Triangulation 

For the algorithm to construct a piecewise affine Lyapunov function we need to 
fix our triangulation. This triangulation is a subdivision of Sj- x W 1 into (n + 1)- 
simplices, such that the intersection of any two different simplices in the subdivision 
is either empty or a /c-simplex, < k < n+ 1, and then its vertices are the common 
vertices of the two different simplices. Such a structure is often referred to as a 
simplicial (n + l)-complex. 

In contrast to |12j . we do not need to have a fine triangulation near 0, but we 
need to ensure that the triangulation respects the periodicity in t of the phase space, 
i.e. it is a triangulation of the cylinder x W 1 . 

For the construction we use the standard orthonormal basis e"o, ex, e~2, ■ ■ ■ , e n of 
M x M n , where e*o denotes the unit vector in i-direction. We also fix a scaling matrix 
S = diag(l, sx, ■ ■ ■ , s n ) with diagonal entries Si > which fixes the ratio of the 
fineness of the triangulation with respect to different directions. We have fixed 
the 1 in t-direction to ensure that the simplicial complex is compatible with the 
periodicity. 

Further, we use the characteristic functions Xj-(i) equal to one if i G J C N 
and equal to zero if i ^ J and the functions R,^ : M n+1 — > W l+1 , defined for every 
Jc{l,2,...,n} by 

n 

R J (x) : =x eo + ^(-l)^ W ^, 
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where x = {xq,x\, . . . ,x n ) and xq = t. Thus Ry (£) puts a minus in front of the 
coordinate X{ of x if % £ J . 

Definition 4.1 Denote by Af the set of all subsets P C x i™ that fulfill: 

i) T> is compact. 

ii) The interior T>° ofD is connected and open. 
Hi) V = V s . 

Note that compactness, connectedness etc. refer to the space x W 1 . 




Figure 1: The triangulation 7~ basic of S^, x R 1 . Note that the operator puts a 
minus-sign in x-directions, but not in the f-direction, which results in the simplices 
shown in this figure. 

Definition 4.2 Let C G N be a given subset of x W 1 . We will define a triangula- 
tion 7~x of a T>k £ N , T>k ^> C, that approximates C. To construct the triangulation 
of a setl~x, we first define the triangulations T basic andT^ slc as intermediate steps. 

1. The triangulation T hasic , cf. Figure [/I consists of the simplices 

&Z,j,a ■= co (ll J + g *«J : J = _1 > 0, 1, 2, . . . , rz^ 

for all z G Nq +1 , all J C {1,2, ... ,n}, and all a 6 S n+ i, where S n+ i denotes 
the set of all permutations of the numbers 0,1,2, ... ,n. 

2. Now scale down the triangulation T bastc with the mapping x, i— > pSx, where 
p := 2~ K T and S is the fixed diagonal matrix defined above. We denote by 
'j-basic ^ g resu m n g se f f ( n _)_ I s ) -simplices, i.e. 

6 :=co(^SR J yz + J2^(i) j = i = -1, 0, 1, 2, . . . ,nj . 
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Note that for each simplex © E T^ asic there is a unique i E No such that 
© E [iT, (i + 1)T] x M n . This follows from the fact that the scaling matrix 
S has 1 as its entry in t-direction and p = 2~ K T . Hence, we can and will 
interpret a simplex T E ~]~^ asic as a set in x M. n in the following step. 

3. As a final step define 

7|:= {ee Tk asic ■ © n C V } 

and set 

V K := (J 6cS}x M n . 
Be?}? 

TVoie i/iaf consists of finitely many simplices due to the fact that C is com- 
pact and the triangulation respects the periodicity. 

Lemma 4.3 Consider the sets C and T>k from the last definition. Then T>k D C 
and T>k E M . 

Proof: T>k is a closed set containing C° and thus contains C E N because C = C° 
by property iii) in Definition |4.1| so C is the smallest closed set containing C°. That 



T>k fulfilles properties i) and iii) of Definition 4.1 follows directly from the fact that 



T>k is a finite union of (n + l)-simplices. To see that property ii) of Definiton 4.1 



is also fulfilled, i.e. that T>° K is connected, notice the following: The definition of 7^ 
implies that for any & E we have &° flC° ^ 0. Because of this, any x E © can 
be connected to a y E C° with a line contained in &° with a possible exception of 
the endpoint x. Because C E N we have that C° is connected so this implies that 



T>° K is connected too, i.e. T>k also fulfilles property ii) of Definition 4.1 and therefore 



T>k E TV. □ 

Remark 4.4 The triangulation T basic is studied in more detail in sections 4-1 and 
4-2 in \22j . A sometimes more intuitive description of &z.j.a is the simplex {x E 
R n+1 : < x CT (o) < • • • < £<j(n) — 1} translated by x i— )■ x + z and then a minus-sign 
is put in front of the i-th entry of the resulting vectors whenever i E J C {1, . . . , n}; 
but no change of sign in the t-coordinate. 

7~K is truly a triangulation, i.e. two different simplices in 7j£ intersect in a com- 



mon face or not at all, as shown in Lemma 4.5 



Lemma 4.5 Consider the set of simplices T£ from Definition 4-2 and let & 



co(xq, xx,..., x n+ i) and &2 = co(yo> Vi, ■ ■ ■ , Vn+i) be two of its simplices. Then 

©i n 6 2 = ©3 := co(z , zx, ■ ■ ■ , z m ), (4.1) 

where zq, z±, . . . , z m are the vertices that are common to ©i and 62, i.e. z% = x a u\ = 
yp(i) f or a, /? E S n+ 2 and i = 0, . . . , m, m E {—1, 0, . . . , n + 1}. 
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Proof: The equation (4.1) follows as in Theorem 4.11 in 



triangulation respects the periodicity. 



and the fact that the 
□ 



One important property of the chosen triangulation is that the simplices are 
sufficiently regular, i.e. the angles all have a lower bound. To make this precise and 
to also measure the influence of K, we prove the following lemma. The matrix Xj^ v , 
as defined in the next lemma, consists of the n + 1 vectors which span the simplex. 
We obtain an estimate on its inverse depending on K. This estimate will later 
be used to estimate the derivative of an affine function on the simplex. 



Lemma 4.6 Using the notation of Definition ^. 4 there is a constant X* , which is 
independent of K and v, such that for all simplices & u E fj^ astc we have 



\X^h < —X* 



2* 
s*T' 



where Xk v 



( (x\ - x ) T \ 

(X2 ~ X ) T 



, xq, . . . , x n+ i are the vertices of & v (in any order) 



\ (x n+1 - x ) T ) 
and s* := min(l, s\, . . . , s n ) > 



Proof: Every simplex in T^ as4C is formed from an (n + l)-simplex 
eo(ico) %)•••) %n+i) ^ 7" basic with corresponding matrix X v . The matrices Xk,v 
and X v relate via 

Xk,u = pX v S. 

Note that for the matrices X u , up to translations, there are finitely many different 
simplices in 7~ basic and also finitely many possibilities of ordering the vertices of any 
such simplex. Hence, there is only a finite number of possibilities of forming such 
a matrix X v . Further, all of them are invertible. This means that we can define 
a 2 > as the minimal eigenvalue of all possible X„ X v . Note that 

^min(S T Xj X^S) > X m i n (X'^ X u ) X m i n (S T S) . 

Hence, 

\\XxlWx < Vn~+l\\X K i\\ 2 



n + l 



< 



Xmin(X^- u XK,u) 

1 n+l 



s*aT 
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as p = 2 K T. This shows the lemma with X* = y/ ^ +1 . Note especially that X* is 
a constant independent of K and v. □ 

In the next lemma we show that a CPA function on a simplicial complex as 
defined above satisfies the assumptions of Lemma 3.3 and the technical assumption 
of Theorem 13. 41 



Lemma 4.7 Let T be a simplicial complex in x W 1 , which is locally finite, 
i.e. each point has a neighborhood U such that U D & ^ only for a finite number 
of simplices © £ T. Denote V = U 6e r©- Let M E C° (£>,§") be a CPA function, 
which is affine on each simplex ofT, and let M(t, x) be positive definite for all (t, x). 
Then M is Lip schitz- continuous on T>° and for each (to,xo) E T>° 

M{{t ,x ) + 6f(to,x ))-M(to,x ) 

hm (4.2) 

0^0+ 9 

exists. 

There also exists (at least) one simplex & u E T and 9* > such that 

(t ,x ) + 9f(t Q ,x Q ) E e u for all 6 E [0,9*]. (4.3) 



M\q u restricted to this simplex is an affine function and the expression in (4-2) is 
equal to 

• M\' & (to,£o) = V;zM|e„(tO) ^o) • f(to,xo), the smooth orbital derivative of the 
affine function M\q v , 

• and M' + (to, xq). 



These expressions are the same for all simplices which satisfy (4.3). 

In particular, M is a Riemannian metric in T>° in the sense of Definition 3.1 
Moreover, 

T 

p(t) T M' + {t,x(t))p{t) dt 

o 

exists and is finite for all solutions x(t) with x(t) E T>° for all t E [0,T] and all 
functions pe C°([0, T], R n ). 



Proof: Let xq := (to,xo) E T>°. Since there are only finitely many simplices 
Si, ... ,6 n E T, which have a non-empty intersection with a neighborhood U of 
xq, and on each of them My is affine and has a finite Lipschitz constant, the overall 
constant can be chosen as the maximum of the finitely many. 

Now we show that there is a 9* > and a simplex & u E T such that xq+9 f(xa) E 
©i, for all 9 E [0, 9*]. Then it is clear by the smooth chain rule that the limit (4.2) 
exists and is equal to the smooth orbital derivative M|g (i, x) = VxM\q v (xo) ■ f(xo) 
of the function M restricted to the simplex & v . Furthermore, (4.2) is equal to 



M' + (xq) by Lemma 3.3 If there are two simplices ©i and ©2 with property (4.3), 
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then M\q 1 (xq) = M\q 2 (xq) and also (M|g 1 )' + (xo) = (M|e 2 )+(^o) by Lemma 3.3 
since M\ 6l (x) = M| 6a (x) for all x = x + 9f(x ) with 9 G [0,9*]. 



Now we show that there exists a simplex with property (4.3). Indeed, there 



is a J G N such that xq + g JJ for all j > J where U is the neighborhood 

of xq from above. Assume that there is no simplex ©„ and no 9* > such that 
xo + Of(xo) e & u for all 9 G [0, 9*}. Then 

|xo + ^M,J> j|n6 fc = 

for all simplices ©i, . . . , ©jv* G T to which xq belongs, since if one such point was 
in &k, then, due to the convexity of &k, the whole line between that point and xo 
would be in & u . 

Since there are infinitely many points in |xo + ^J \ J — j\ ; which are in U, but 
only finitely many simplices that have nonempty intersection with U by assumption, 
at least one of them, say & v , must contain infinitely many such points. Since the 
sequence xq + ^y 1 ^ converges to xq as j — )■ oo and <5 V is closed, xq G & u and hence 
Sjy must be one of the &i, . . . ,&n* defined above, which is a contradiction. This 
shows the statement. 

For the last statement, we show that the function 9 — > M^_(SqXq) is RHS contin- 
uous. Then the function p(t) T M' + (t, x(t))p(i) is RHS continuous and bounded, and 
thus it is integrable. 

Fix xo and a neighborhood U such that there are only finitely many simplices 
with nonempty intersection with U. Denote by &i, . . . ,&n the subset of those 
finitely many simplices such that for each ©j there is a 9* with xq + 9/(xq) G &i for 
all 9 G [0,0*], these are the ones the qualify to be &{xq). 

Now take a sequence 9k — > + and we seek to prove that M+(Sg k xo) — > M' + {xo). 
Since Sg k xo G U, if k is large enough, there are only finitely many simplices which 
contain infinitely many elements Sg k xo of the sequence. We show that these simplices 



are in fact a subset of ©i, . . . , ©jv as defined above, satisfying property (4.3). If this 
was not true then there would be a simplex & and a sequence (the subsequence 
from above) of points Sg k XQ G 6, xq G ©, but xq + 9f(xo) © for all 9 > 0. The 
simplex © is the intersection of finitely many half-spaces. The point xq lies on the 
boundary of (at least) one of these half-spaces, since otherwise a neighborhood of 
xq would lie in 6. Hence, there is a half-space of the form {x \ (x — xq) ■ h < 0} 
with n G M n+1 \ {0} such that /(xo) ■ n> 0. Hence there is an e such that v ■ h > 
for all v with \\v — /(xo)|| < e - Now we have 

© 3 Sg k x = x + 9 k / f(x + \(Se k x - xq)) dX. 
Jo 

If k is large enough, then v := J Q /(xo + X(Sg k xo — xq)) dX fulfills — /(xo)|| < e by 
continuity of / and thus Sg k XQ G" 6, which is a contradiction. This shows that each 
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simplex which contains infinitely many elements of the sequence fulfills property 
(4.3) and thus is a candidate for <&{xq). 

Let 6 be one of the simplices which contains infinitely many elements of the 
sequence and define the subsequence 6k t by choosing the next 9k with the property 
Se k xo G S. We do the same for all (finitely many) simplices that contain infinitely 
many elements. For the convergence of the overall sequence it is enough to show that 
every subsequence Sg k XQ converges to the same limit, since there are only finitely 
many. 

For each subsequence, the elements are in one simplex and on this simplex M' + 
is a continuous function, so it converges. The limit is the same, as M' + {xq) is the 
same no matter which simplex &(xq) we choose. □ 



Remark 4.8 The condition locally finite for the triangulation is indispensable as 
shown by the following example: Let /(l/2,0) = (0, 1) and for every n G N define 
the triangle & n := co((0, 0), (1, l/n), (1, l/(n + 1))). Then clearly there is no n G N 
with a corresponding 9* > such that (1/2,0) + 5/(1/2,0) = (1/2, (9) G & n for all 
9 G [0,9*]. 



4.2 The semidefmite optimization problem 

For each simplex & u G 7^ we denote 

h u : = diam(©„) = max \\x — y\\2- 

Note that for our triangulation we have with (5„ = co(xo, ■ ■ • , x n +i) and S* 
y/n + 1 max(l, s%, . . . , s n ) the estimate 



h v = max \\xk — X1W2 < S* p = S*2 T. 
fc,ze{o,...,n+i} 



Moreover, denote 



max 
x£G v ,i,j£{0,...,n} 



max 

a:ES„,i,j,fc£{0,...,rt} 



d 2 f(, 



dxidxj 



d 3 f(j 



dxidxjdxk 



where xq := t 



if / G C 3 



(4.4) 

(4.5) 
(4.6) 



Variables 

The variables of the semidefinite optimization problem are 

1. Mjj(xfc) G K for all 1 < i < j < n and all vertices x^ of all simplices G v = 
co(xo, • • • , x n+ i) G - values of the Riemannian metric at vertices 

2. C v G M.q for all simplices & u G - bound on M in 6^ 

3. D u G Mjj" for all simplices 6 F G 7^ - bound on derivative of M in 6 U 
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Thus we have 2s + ^n(n + l)v variables, where s denotes the number of simplices 
and v the number of vertices. 

In this section, there is no objective function, so we are considering a feasibility 
problem - see Section [4. 4| for a suitable objective function. The equality Constraint 
1. can be incorporated by choosing the variables in this way and thus is no actual 
constraint. Note that Constraint 3. is linear and Constraints 2., 4. and 5. are 
semidefinite. Constraints 2. and 5. need to be satisfied for each simplex and 
then for each vertex, so vertices common to more simplices need to satisfy several 
constraints. In Constraint 4., however, each vertex only needs to be checked once. 
Note, however, that we can replace the individual constants C v for each simplex by 
their maximum C, so that also Constraint 2. only needs to be checked at each vertex. 
Note that if the triangulation is fine enough and the system has an exponentially 
stable periodic orbit, this more restrictive form of the constraint can be satisfied, cf. 



the proof of Theorem 5.1 
1. Periodicity 



M ij (0,x k ) = M ij (T,x k ) 

for all 1 < i < j < n and for all vertices at times and T, i.e. for all vertices 
(0,Xk) and (T,x k )- Note that by construction of the triangulation (0,x k ) is a 
vertex of a simplex in 7]£ if and only if (T, x k ) is. 

2. Bound on M 



M(x k ) * CJ 

for all vertices x k of all simplices & u = co(xq, . . . , a? n +i) G 7]£, where the 
symmetric matrix M{x k ) is defined by setting Mji(x k ) '■= Mij(x k ) for all 
1 < i < j < n. 

3. Bound on derivative of M 



J n+l 

for all I = 0, . . . , n, 1 < i < j < n and for all simplices & u E TS, where 



W; 



VjMjjlg (x) for all x G & v = co(xq, • • • , x n +i), which is given by 



/ Mijixx) -Mij{xo) 



w 



'j 



X k]v 



m+l 



(4.7) 



\ Mij(x n +l) ~ Mij(xo) 



where Xk v 



( (x x - x ) T \ 

(x 2 ~ X ) T 

\ (x n+1 - x ) T J 



9 (n+l) x (n+l) 
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4. Positive definiteness of M 

We fix e > 0. 

M(x k ) h e I 
for all vertices x k G of all simplices & v G 7j£. 

5. Contraction of the metric 



M(x k )D x f(x k ) + D x /(x fc ) T M(x fc ) + «j • f(x k ))ij=i,..., n + (E„ + 1)1 ± 

for all simplices ©„ = co(xo, . . . , £En+l) G 7]£ and all of its vertices k = 

0, ...,n + 1. Here, f(x) = f ^ and (u;?. • /(5fc))i,j=i,...,n denotes the 

symmetric (n x n) matrix with entries ■ f(x k ), where was defined in 
(4.7) and is the same vector for all vertices in one simplex and 

f hunB v yn + \huD v + 2n(n + 1)C„] if / G C 2 

\ hXv^+Ttl + 4^5^ + 2n(n + l)B^ v C v ] if / G C 3 1 j 

Remark 4.9 In Constraint 3. we have claimed that the gradient of the affine func- 
tion Mij\ & , i.e. VxMij\ & = w"j, is given by the expression in [4.1) ■ For a proof of 
this fact and, moreover, that the definition is independent of the choice of the vertex 
xq, see fT2\ Remark 2.9]. 

Remark 4.10 The constraints above are easily transferred into the standard form 
Y^i=\FiDi — Fq >z in the following way: Denote by y±, . . . , y n (n+i)/2 the matrix 
elements Mij(x\), 1 < i < j < n, by the following elements yi the matrix 

elements of Mijijb?) etc. and finish the vector y by the C v and D u . This results in 

n i n(n+l) i 

m = Is H — — -v as above. 

Now Constraints 2. and 4- are expressed in (n + 2)s, v blocks of size n each, 
respectively, in the matrices F^. Constraint 3. is expressed in 2 ■ "( n 2 +1 ) conditions 
(for each i, j) for each of the n + 1 entries of the vector w\y This needs to be 
considered for each simplex, so that we have n(n + l) 2 s blocks of size 1. 

Finally, Constraint 5. is expressed in (n + 2)s blocks of size n, since for every 
simplex each of its vertices needs to be considered. Note that E v depends linearly on 
C v and D u and depends linearly on Mij(xk). 

The size of the matrices Fj is thus a block- diagonal structure with n{n + l) 2 s 
blocks of size 1 and v + 2(n + 2)s blocks of size n. 

Remark 4.11 Note that Constraint 2. implies that niaXfc = o n +i \Mn(xk)\ < 
1 1 M(x~k) 1 1 max < ||-^(5fe)||2 < C u since M(x k ) is positive definite. Note also that 
Constraint 3. implies that < D v . Moreover, Constraint 5. is equivalent to 

X max (M(x h )D x f(x k ) + D x f(i k ) T M{x k ) + (w% ■ f{x k ))ij= x ,...^ + E v < -1 
where X ma x denotes the maximal eigenvalue. 
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4.3 Feasible solution is CPA contraction metric 



A solution of the semidefinite optimization problem returns a matrix Mij(x k ) at 
each vertex x k of the triangulation for 1 < i < j < n. We define the CPA metric by 
affine interpolation on each simplex. 

Definition 4.12 FixC and a triangulation withT>K = Ue„eT c ®«" Let Mij(x k ) 
be defined by a feasible solution of the semidefinite optimization problem. Let (t,x) = 
x G & u = co(xo, • • • , x n +i) such that x = Y^k=o ^k%k with Afc G [0, 1] and Ylk=o 
1 . Then define 

n+l 



M(x) =^2x k M(x k ) 



k=0 



Lemma 4.13 The matrix M(x) as in Definition is symmetric and positive 



definite for all x G T>k- The function M(t,x) is periodic in t with period T. 



Proof: The symmetry follows directly from the symmetry of M(x k ): 

n+l n+l 

Mij[x) =J2 X k M ij{xk) = ^2X k M ji {x k ) = Mji(x). 
k=0 k=0 

The positive definiteness also follows from the positive definiteness of M(x k ), 
using that the minimal eigenvalue A m i n is a concave function. Indeed, consider 
x G & u = co(xo) • • • Note that X m i n (M(x k )) > eg for all k = 0, . . . ,n + 1 due 

to Constraint 4. Thus, 

/n+l \ 

X min (M(x)) = \ min AfcM(x fc ) 

\k=0 J 

n+l 

> ^AfcA mfa (M(x fc )) 

k=0 
n+l 

> eo A fc 

k=0 

= eo>0. 

The T-periodicity follows directly from the definition and the triangulation: 

n+l n+l 

M (0, x) = J2 X kM{0, x k ) = A fc M(T, x k ) = M(T, x). 

k=0 k=0 

□ 

We will now relate M'(x) to M'(x k ), as well as M(x)D x f(x) to M(x k )D x f(x k ). 
For the proof we will need the following auxiliary result, see [21 Proposition 4.1 and 
Corollary 4.3]. 
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Lemma 4.14 Let f G C 2 (R n+1 ,W l ) and x G co(x , . . . , x n +l) = & u with x 
Y2to AfcX fc , A fc G [0, 1] and ^k = 1 - 

Then for all I G {1, . . . , n} we have 



n+l 



especially 



k=0 



n+l 



< max \\H fl (x)\\ 2 hl, 



Xk) 



k=0 



where h u = diam((S t ,), B u = maxg ee!/j j )je r 0) n i 



, xq := * ari( ^ Hf^x) :- 



( d 2 h(x) 
y dxidxj 



3 ' i J=0,...,n 



denotes the Hessian of ft. 



In the following we restrict ourselves to one simplex & u G Tj£. First we need to 
define M'(x) (the orbital derivative) for a general point in the simplex & u . Note 
that for all points of a simplex the vector V^M^-fx) is the same, but the contribution 
of f{x) is different. 



Definition 4.15 Let Mix) be as in Definition 4-12, Fix a point x G T>° K and a 
simplex G v = eo(xo, • • ■ , %n+l) £ 7k suc -h that x + 9 fix) G & v for all 9 G [0, 9*], cf. 
Lemma ^.7 Then x = J2k=o ^k x k with X k G [0, 1], X^fcio ^k = 1 an d 

(M,, )'+(*) = M'^ix) = V £ M tj \ 6ij ix) ■ fix) = w% ■ fix). 

Note that My| 6 is an affine function and its gradient was defined in Constraint 
3. and is the same vector for all points x G & u . Note that for vertices x = Xk this 
term appears in Constraint 5. 



Lemma 4.16 Let f G C 2 and let Mix) be as in Definition 4-12 Fix a point x G T>° K 
and a corresponding simplex G u G T K as in Definition 4-15 
Then we have the following estimates for all x G & u 



n+l 



iM ij y + ix)-Y,\< j -h^k) 

k=0 

n+l 

{M{x)D x f{x))ij -J2^iMix k )D x fixk))ij 

k=0 
n+l 

iD x fix) T Mix)) i3 -Y, X kiD x fixk) T Mix k )) ij 



k=0 



< Vn + lB„D v hl 

< n(n + \)B v C v h v 

< nin + l)B v C u h v , 



where B u = niax i66 „,i J( :{o,.,n} 



d 2 f(x) 



, xq := t and h v = diam((3y). 
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Altogether we have 



(M(x)D x f(x) + D x f{x) T M{i) + M' + (x)) 

~ E£o X k (M(x k )D x f(x k ) + D x f(x k ) T M{x k ) + (<• • /(x fc )) 4 , i= i,. 



< h u B u (Vn + Dh,A, + 2n(n + l)C v ) = — 

n 

If f G C 3 we obtain in addition the estimates 

n+l 

(M(x)D x f(x))ij -J2MM(x k )D x f(x k )) 

k=0 

< n(2jn~+\B v D v + (n + \)B^ v C v )h 



n+l 



(D x ./(5) T M(5)) lj - ^ \ k {D x f(x k ) T M(x k ))i 



k=0 



< n(2Vn + \B V D V + (n + \)B^ v C v )h 



where B 3iV = max ig6j/)ij>e{0v .. ]n} 



d 3 /(*) 



, xq := i. Altogether we have 



(M(x)D x f(x) + D x f(x) T M(S:) + M' + (x)) 

~ X k (M (x*)^/^) + D x f(x k ) T M(x k ) + (<• • f(x k )) iJ=1> . 



< hl[Vn~Tl{l + 4n)B u D u + 2n(n + l)B 3 , u C u } 



n 



Proof: Step 1: M' 

Fix one simplex (3^ G of the triangulation and let x G <5„. By defini- 
tion of Mij(x) as a CPA function, interpolating My(%) at the vertices, we have 
V x Mij(x) = w%j for all x G <5j,, where M is restricted to the simplex, cf. Definition 
Thus, letting x = Efcio X k%k with \ k G [0, 1] and Y^k=o X k = 1 



4.15 



V*M«L (5)-/(x) 



'n+l \ / n+l N 

x khxk) + <■ • /(x) - £ x kf(*k) 

\k=0 J \ k=0 / 
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Hence 



n+1 



(My);(x)-^A fc ^7(z fe ) 



k=0 



since 
is 0. 



fix) ~ Efcio x kf(x k )) 



< IK- Hi 



- E5S x kf(x k )) 



n+1 



/(^) - X kf{xk) 



k=0 



as the t-component 



Now we use Lemma 4.14 for the C 2 function /, establishing that ||/(x) — 
ZSWO^Iloo < {n + l)B u h 2 v . Thus, 



n+1 



{M ij )' + {x)-Y J x kW v ij -f{x k ) 



k=0 



< ||<||i(n + l)B v hl 



Using that ||u&||i < D u holds from Remark 



4.11 



131 

we obtain 



n+1 



(M ij )' + (x)-Y, x kW^ j -f(x k ) 



k=0 



< {n + l)D v B v h 2 v . 



Step 2: MD x f 

We consider (M(x)D x f(x))ij = J2?=i Mu(x)(D x f(x))ij. We first consider two 
scalar-valued functions g and h, where g{x) = J2k=o x d( x k) an d h is C 1 in x. We 
have, using Taylor expansion for h at x^, i.e. h(x) = h(x k ) + V/i(x*)(x — x k ), where 
x* lies on the straight line between x k and x, 



n+1 



g(x)h(x) 



^2 X k9(x k )Kx) 

k=0 

n+1 

}. ] Afcff(gfe)[fc(x fc ) + V^/i(x*)(x - x fc )] 



fc=0 



n+1 



g(x)h(x) - ^2 x kg(xk)h(x k ) 



k=0 



< max ||V £ /i(x*)||i/ij, max |#(xfc)| 
z*e6„ fe=o,...,n+i 



where ||x — x k \\oo < \\x — x k \\2 < /i^. 

Applying this to g(x) = Mu(x) and /i(x) = (D x f(x))ij we obtain with 
max fc=0 ,...,n+i |Mij(£jk)| < C v from Remark |4~TT| and || V^(D x f(x*))ij\\i < (n + 



1) max^gg^ j^gjo jn j )fc ^ 



< (n + !)-£?!, for all Z, j G {1, . . . , n}, using that 



x, Xfc G and © u is a convex set. Thus, 

1 n n+1 



/=1 fc=0 



< n(n + \)B v h v C v 
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Hence, 



n+l 



(M(x)D x f(x))ij -J2 X k( M (xk)D x f(x k ))i 



k=o 



< n(n + \)B v h v C v 



A similar estimate holds for D x f(x) T M(x). 
Step 2': MD x f 

If / G C 3 , then we can derive an estimate for this term which establishes order 
h 2 . We consider two scalar-valued functions g and h, whe re gi x) = J2k=o ^kg(x k ) 
for x = X^fcio ^k%k and /i(ir) is C 2 in i. We apply Lemma 4.14 to (g • h), yielding 



n+l 



g{x)h(x) - ^2 hg(x k )h(x k ) 



k=0 



< max\\H(x)\\ 2 h* 



where the matrix H(x) is defined by (H(x)) 



km 



^h + gJ^- and 



d 2 (g-h)(x) 
dx k dx m 



Note that -J- 

OX r , 



(4.9) 

-(g-h) = 



2 



dx k dx r 



<g-h) 



d 2 g , dg dh dg dh d 2 h 
h + 7; — 7; h 7; — 7; h 5- 



dx k dx 



dx m dx k dx k dx r 



dxudxr 



Applying this to g(x) = Mu(x), we observe that, since g(x) = Mu(x) is affine on 



dg 



the simplex, ^ 
d 2 



w 



ilJm 





and Q Xk Q Xm (x) = for all x G & u . Hence, 



dxfcdxj 



, , w~n , dh(x) . ... dh(x) _ _ ._. d 2 h(x) 
(g ■ h)(x) = ( w ») m ^-l + K)*-^ + M«(x^ 



Using h(x) = (D x f(x))ij, we obtain with ^ - g^^- 



and 



dx k dx m ' 
d 2 h _ d 3 h 



dxkdx m dxj ' 



j ^ 0, defining 5 3ji , = max 5eei/ji)J y 6{0j ... jn} 

d\g-h){x) 



9 :i f(x) 



dxidxjdxi 



\{H(x)) km \ 



dx k dx m 

< \(v%) m \B v + \(wZ) k \B v + B 3tU \M a (x) 



Thus, using \\H\ + H2II2 < ll-^ilb + ll^lb, as well as ||#i||2 < V n + l||i/i||i, 
||-ff 2 ||2 < v^Tlll^lloo and ||i?|| 2 < (n + max we obtain 



\H\\ 2 < 2,/n + l\\v%\\ 1 B v + (n + 1)B 3; y max max | M# (x) \ 

:r£6„ l<i<l<n 

< 2V^+lD u B u + (n + l)B 3iI/ C v , 



using Remark 4.11 . Hence, (4.9) establishes 



n n+l 



Z=l i=l fc=0 

< nh 2 u {2^/^+lD v B v + (n + l^C). 
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which proves the lemma. □ 

Now we can estimate the value of Lm for all points x for the CPA metric M, 
given our constraints on the vertices. 

Lemma 4.17 Let all constraints be satisfied. Then the CPA metric M defined in 



Definition fulfills: 

\ m ax(M(x)D x f(x) + D x f(xfM(x)+M^(x)) < -1 (4.10) 
for all x G V° K . 

Proof: The maximal eigenvalue is a convex and thus sublinear function, i.e. for 
L, S G S n symmetric (but not necessarily positive definite), we have 

^max{L + S) < X ma x (L) + X max (S). (4.11) 

We show that \\ m ax(S)\ < \\S\\2 holds for a symmetric (but not necessarily 
positive definite) matrix S G M nxn . Indeed, denote by Ai < . . . < X n the eigenvalues 
of the symmetric matrix S with corresponding eigenvectors vi, . . . ,v n , forming a 
basis of R n . Then 

Svi = XiVi =^ S T Svi = SXiVi = Afuj. 
Thus, Vi is an eigenvector of S T S with eigenvalue A^. Hence, 



= max(|Ai|, . . . , |A„|) 

= max(|Ai|, |A n |) > |A n | = \X m ax(S)\. (4.12) 
Fix x e V\ and 6„ £ 7| such that x + 6f(x) G & u for all 9 G [0, 6*} as in 



Definition 



4.15 



Then x = YJk=o Afc^fc, YJk=Q Afc = 1 and X k G [0, 1]. Now, using 
(I4TTT1) and ||5|| 2 < max, where HSHmaa; = maxi< i<j< n \ Sij\ for a matrix S G S n , 
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we have with Remark 14.111 



\ max {M(x)D x f(x) + D x f(xY M(x) + M' + (x)) 

/n+l 

< \max Yl Xk (M(x k )D x f(x k ) + D x f{ Xk ) T M(x k ) + (v% • f(£ k ) )< j=l,..., n 



\k=0 



+X max ( M{x)D x f{x) + D x f(x) T M(x) + M' + (x) 
n+l 

- Afc (M(x k )D x f{x k ) + D x f{x k ) T M(x k ) + • /(%)) 

fc=0 

< ^ [M(x k )D x f(x k ) + D x f{x k ) T M{x k ) + (w^ ■ f{x k ))ij=i,..., n 



i,j=l,...,n 



n+l 



fc=0 



+ 



M{x)D x f{x) + Af(x) + M + (x) 

n+l 



k=0 



- ^ A fc (M(x k )D x f(x k ) + D x f(x k ) T M(x k ) + (w£ • /(%))i,i=i, 



n+l 



< ^A fc (-l-^) + 



fc=0 



M(x)D x f(x) + D x f(x) T M(x) + Mi(x) 



n+l 



- A fc (M(x k )D x f(x k ) + D x f(x k ) T M(x k ) + • /(x fc )) M =i,...,n 

fc=0 

< -1 - ^ + ^ 
= -1 



by Lemma 4.16 and the definition of E v . 

We summarize the results of this section in the following theorem. 



□ 



Theorem 4.18 Let be a triangulation as in Definition and T>k '■= 
Ue„eT c Let all constraints be satisfied. Then the CPA metric defined in Defi- 
nition \4-l£\ fulfills: 

1. M(x) is symmetric, positive definite and a T -periodic function and thus defines 



a Riemannian metric on T>° K in the sense of Definition 3.1 



2. Lm(x) < — 2^T — — ~W ^ for all x G D° K , where fj, T , 



max x( zD K \ max {M(x)) > 0, Lm is defined in Theorem 3.4 and C :- 
3. If G C T>° K is positively invariant, then 



p{t) T M' + (t,x{t))p(t) dt 



24 



exists and is finite for all solutions x(t) with x(0) G G and all functions p G 
C°([0,T],R n ). 



Hence, M satisfies all assumptions of Theorem 3.4 



Proof: Part 1. follows directly from Lemma 4,13 and 4.7 



For Part 2., note that [i max exists and is positive, since M is positive definite for 
all x and depends continuously on x. Then 



L M (x) 



w 



sup 

weK n \{o} 



[M{x)D x f{x) + \M' + 



X )\W 



< 



< 



w T M(x)w 

1 \mg X (M(x)D x f(x) + D x f{x) T M(x) + M' + (x)) 

2 Xmax{M(x)) 



using Lemma 4.17. Moreover, Constraint 2. yields that \i max < max e G7 -c C v < C. 

Part 3. follows from Lemma 4.7 using that G C T>° K is positively invariant. This 
shows the theorem. □ 



4.4 Objective function 

While we are primarily interested in the calculation of any Riemannian contraction 
metric, i.e. a feasible solution of the semidefinite optimization problem, the optimiza- 



tion problem also allows for an objective function. Theorem 4.18 suggests a possible 
objective function, namely simply the maximum over all C u , i.e. C = max 6 ^ g j-c C v . 
We can either implement this by adding C as an additional variable with constraints 

C V <C 

for all v. Or we can replace the variables C u in Constraint 2. by the uniform bound 
C directly, which reduces the number of Constraints 2. to v, the number of vertices. 
By minimizing the constant C, we minimize the bound on the maximal Floquet 
exponent. Note, however, that C has a lower bound given by eo- This means that 
by choosing eo too large, the estimate on the Floquet exponent will be very rough. 



5 Feasibility of the semidefinite optimization problem 

In the next theorem we assume that there exists an exponentially stable periodic 
orbit. Then we can show that the semidefinite optimization problem has a feasible 
solution and we can thus construct a suitable Riemannian metric. We have to 
assume that the triangulation is fine enough, i.e. in practice we start with a coarse 
triangulation and refine until we obtain a solution. The triangulation has to stay 
suitably regular, i.e. the angles in simplices have lower and upper bounds. For 
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simplicity, we use the reference simplicial complex and scale it uniformly as described 



in Section 4.1 but other refinements are also possible. 



Theorem 5.1 Let the system x = f(t,x), f G C 2 (S^ 



have an exponen- 



tially stable periodic orbit Q with basin of attraction A(Q). Let C C A(£l), C G N , 



$7 C C be a compact set in the cylinder 



Fix eo > 0. 



Then there is a K* G N such that the semidefinite optimization problem is feasible 



for all triangulations 7j£ as described in Section 4-1 
choose the constants C u and D 
D u = D for all v. 



with K > K* . Note that we can 
to be the same for each simplex, i.e. C v = C and 



Proof: 

1. Smooth Riemannian metric 

Denote the maximal real part of the Floquet exponents of the exponentially 
stable periodic orbit Q by — vq < and set e := ^. 

Since C C A(Q) is compact and A(£l) open, there is a positive Euclidean 
distance between C and the boundary of A(Q). Let d > denote this distance 
if it is finite and otherwise set d := 1. Now define C* to be the set of all 
W 1 that have Euclidean distance less than or equal to d/2 to C. 
Clearly C* £ N . For all large enough K 6 No the Euclidean distance from the 
boundary of T>k '■= &u to C is bounded above by max 6 ^ g7 -c diam(6^) < 



S*2- K T, see (|44|). Thus, there is a K** G N such that C C V K C C* for 
every K > K**. 

Now apply the following theorem to the above defined C* and e. 

Theorem 5.2 (Theorem 4.2, [8j) Assume that f G C°{R x IR^IR™) is a 
periodic function in t with period T, and all partial derivatives of order one with 
respect to x are continuous functions of (t,x). Let £1 := {(t,x(t)) 6 Sj, x IR n } 
be an exponentially asymptotically stable periodic orbit, A(Q) be its basin of 
attraction, and let the maximal real part of the Floquet exponents be —vq < 0. 
Then for all e > and all compact sets C* with £1 C C* C A(£l) there exists 
a Riemannian metric M £ C 1 (C* ,M. n><n ) , such that Lj^(t,x) < —uq + e holds 
for all (t,x) GC*. 

Remark 5.3 The proof of this theorem, cf. J21 Theorem shows that we 
have M G C 2 {C* ,R nxn ) , if f G C 2 {S^ x M n ,M n ) ; as is the case by our as- 
sumptions. 

Since both M(x) and \ m in are continuous functions and M is positive definite 
for all x, G C*, there exists an e\ > such that 

X min (M(x)) > ei 
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for all x G C* , as C* is compact. 

Moreover, since Lj^ is a continuous function (note that M is smooth) satisfying 
Lj^{x) < — y for all i £ C*, we have 

\ m a X (M(x)D x f(x)+D x f(x) T M{x) + M'(x)) < 

for all x G C*. As M(x)D x f(x) + D x f{x) T M{x) + M'(x) and A max are con- 
tinuous functions, there exists an £2 > such that 

Xmax(M(x)D x f{x) + D x f(x) T M(x) + M'(x)) < -e 2 

for all x G C*, as C* is compact. 

Now define M(x) := max (f^, |0 M(x). Then 



A min (M(x)) > e 



(5.1) 
(5.2) 



A maa; (M(x) J D :l ./(x) + Z) x ./(x) T M(x) + M , (x)) < -2 
for all x G C* . 

M is C 2 on the compact set C*, so that we can define the following constants 



M * 



M 2 * 



max\ max (M(x)) 

xGC* 

max max ||VgMy(:r)||i 

x£C* l<i<j<n 

max max H-ff^^x)]^ 

xeC* l<i<j<n J 



where i?My(^) denotes the Hessian of Mij(x). 

2. Assigning the variables of the optimization problem 

Set 



F : = max||/(x)|| c 

zee* 



B : = max 

xeC*,i,je{o,...,n} 



£?3 := max 

sec* ,i,j,fee{0,...,n} 



d 2 /(x) 



dxidxj 



dxidxjdxk 



if / G C 3 



Finally, we define, using the constant A* from Lemma 4.6 

' S*X* 



C* 

c 

D 

K 

hi 



M 2 *(re + 1) 



2s* 



+ Vn+1 



M* 

(C* + M*)(n + l) 

[2nB(s/n + lD + 2n{n + l)C\- 1 if / G C 2 

[2n( v / ^TT(l + An)BD + 2n(n + l)B-zC)}- 1 / 2 if / G C 3 
1 

2nFC* ' 



27 



Now define 



K* := max 



ln(5*r)-ln(min(^,/i*,l)) 



In 2 



,K* 



let K > K* and consider the triangulation T^. Note that, since K > K* , we 
have by (Ol) 



h v < S*2- K T < S*2~ K *T < mm(hl, h* 2 , 1). 



(5.3) 



Interpolate M(x) on this triangulation, i.e. assign the variables Mij(x k ) for 
all vertices x k of the triangulation; moreover, set C u = C and D u = D for all 
v. By (5.3) and the definition of h\ this ensures 

1 



E v < 



(5.4) 



3. Auxiliary results 

To check the feasibility in the next step, we need some estimates. 

Co nside r a simplex & v = co(xo, xi, ■ ■ ■ , x n +i) G 7]£. Denote, as in Constraint 

/ (xi - x ) T \ 
(x 2 - x ) T 

3, \4.7\ w% := X£ v M* ju , where X KtV = . 

\ (X n +1 ~ X ) T J 



and 



M, 



/ Mijixx) -Mijixo) \ 
Mij (x 2 ) -Mij (x ) 

\Mij(x n+1 )-Mij(xo)J 



(5.5) 



Note that M{j(x) is two times continuously differentiable on & u C C* and for 
i, j G {1, 2, . . . ,n + 1} we thus have by Taylor's theorem 

Mij(x k ) = Mij(x ) + ViMy(xo) • (x fc - x ) + ^(x fc - 5 ) r H Mij (z k )(x k - x ), 

where Hm^ is the Hessian of My and % = xo+$fc(xfc — xo) for some G ]0, 1[. 
By rearranging terms and combining, this delivers with (5.5) 



M^-XK^VzMijixo 



We have 



/ (xi - Xo) T H Mij (zi)(xi - x ) \ 
(x 2 - x ) T H Mi:i (z2)(x2 - Xq) 

\(x n +l - XQ) T H Mi ,(z n+ \)(x n+ l - X )J 



(5.6) 



|(x fc - xq) 1 H Mij (zk)(x k - x )\ < lit\\H Mi:j (z k )h < M%hl. 
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Note that z k G & u since the simplex is convex. Hence, by (5.6), 



\M, 



(5.7) 



Now we need to obtain an estimate on V xMij{x k ). For k G {1, 2, . . . , n + 1}, 
/ G {0, . . . , n}, there is a z k i on the line segment between x k and xo, such that 

diMij(xk) - diMij(xo) = VxdiMij(z kl ) ■ (x k - x ), 

where d\M^ denotes the l-ih component of V^Mjj and d[Mij is the gradient 
of this function. Then we have 



\diMij(x k ) -diMij{x )\ < \\VxdiMij(z k i)h\\xk - x \\ 2 < Vn + lM%h v 
Hence, 

\\VzMyixk)- ViMij{x )\\x < (n + 1) 3/2 M 2 *V 



From this, Lemma 4.6 and (5.7) we obtain the estimate 



K--V 5 Afi i (5 fc )||i 

< \\X K ] u M* jtU - VxM^xo)^ + \\VzMijixo) - VzMijixk)^ 

< WX^UM*^ - Xk^VsM^xo)^ + (n + l) 3/2 M 2 */i„ 
2 X X* n + 1 



s*T 2 
S*** n + 1 

~s* 2~~ 



M 2 */i 2 + (n + l) 3/2 M 2 *^ 
Af|ft„ + (n + l) 3 / 2 M£/i„ 



(5.8) 



using (4.4) and the definition of C*. 



A further useful consequence, which we need later, is that 

< C*h v + M{. (5.9) 

4. Feasibility 

Now we check that the assignment in Step 2. is a feasible point of the semidef- 
inite optimization problem. We discuss each constraint in the following. 

Constraint 1. This follows directly from the periodicity of M. 
Constraint 2. This follows directly from the definition of C. 
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Constraint 3. We have for all 1 < i < j < n and any I = 0, . . . , n by ( 5.9 ) 



IKOil < IKIli = Wx^m^iu < c*k + m* < 

by definition of D as h v < 1. 



n + 1 



Constraint 4. This follows from (5.1). 

Constraint 5. This is the main step and will be shown below. 

We fix a simplex 6 U = co(xq, • • • , x n+ \) E and a vertex x^, k E {0, . . . , n + 



1}. With Remark 4.11 we need to show that 
Xmax (M(x k )D x f(x k ) + D x f{x k ) T M{x k ) + «, • f(x k ))u =1 



j=l,...,n 



+ E U < -1. 



We have E v < cf. (5.4). Hence, using the sublinearity of A max , we have 



Xmax (M{x k )D x f(x k ) + D x f{x k y M(x k ) + (wfj ■ f(x k )) i>j= i ; ... >n j + E v 

< Xmax {M{x k )D x f(x k ) + D x f{x k ) T M{x k ) + M'(x k )) 

+Xmax{{Wij ■ f{xk))i,j=l,...,n ~ M' (x k )) + - 

< -2 + - + n max \v% ■ f(x k ) - Mlj(x k )\ 



where we have used (5.2) and Xmax{S) < ||51|2 5; ^ 1 1 S \ \ max : 

cf. RVM. 

Using M^(x fc ) = V x Mij(x k ) ■ /(%), we obtain 



A ma * (M^Ac/C^) + D x f{x k y M{x k ) + • /(x*))i,j=l,...,n) + 
3 



< 
< 
< 



+ n max ||u£- - VzMij(£ fc )||i • ||/(%) 

l<i<j<n 

+ nC*h v F 



2 

3 1 
2 + 2 



-1 



with (5.8), the definition of F and using h v <h\ = 2n FC* • 

This proves that the constraints are fulfilled and the optimization problem has a 
feasible solution. □ 

Acknowledgement The first author would like to thank Michael Hintermiiller and 
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